THE GENERALIZED MADER'S INVERSION 
FORMULAS FOR THE RADON TRANSFORMS 



Y. A. ANTIPOV AND B. RUBIN 

Abstract. In 1927 Philomena Mader derived elegant inversion 
formulas for the hyperplane Radon transform on R™. These formu- 
las differ from the original ones by Radon and seem to be forgotten. 
We generalize Mader's formulas to totally geodesic Radon trans- 
forms in any dimension on arbitrary constant curvature space. An- 
other new interesting inversion formula for the A;-plane transform 
was presented in the recent book "Integral geometry and Radon 
transform" by S. Helgason. We extend this formula to arbitrary 
constant curvature space. The paper combines tools of integral 
geometry and complex analysis. 



1. Introduction 

Let g{9,s) = J f(x) dm(x) be the Radon transform of a sufficiently 
good function / on the Euclidean space 1R™. Here, the integral is taken 
over a hyperplane x ■ 9 = s, where x ■ 9 = x\9\ + . . . + x n 9 n is the 
usual inner product, 9 is a point of the unit sphere S n_1 = {x 6 I" : 
x\ + . . . + x 2 n = 1}, and dm(x) stands for the volume element on the 
hyperplane x-9 = s. Mader's inversion formulas [M] have the following 



elegant form: 



f(x) = A — F (x,t)\ t=0 if wis even, (1.1) 
f(x) = A x — F 1 (x,t)\ t=0 if n is odd, (1.2) 



where 



(_l)("-2)/2 (_l)("-l)/2 

A) = — ; rr; , A\ — 



n{n-2)\a n - 2 2(n-2)!a„_ 2 
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Fo(x,t)= / G(x, s) log I s—t I ds, F 1 (x,t)= / G(x, s) sgn(s — t) ds 



G(x,s) = / g(6,s + x-9)da(6), 

gn-l 

do~(6) is the surface element, and cr n -i is the surface area of E> n ~ 1 . 

An interesting feature of Mader's result is that she used tools of 
complex analysis (unlike most of other authors in the area). 

In the present paper we derive more general inversion formulas, hav- 
ing the same nice structure as (11 .ip and (ll.2p and applicable to totally 
geodesic Radon transforms in any dimension on arbitrary constant cur- 
vature space. We also generalize one recent inversion formula by S. 
Helgason, which has a close flavor and is described below. 

Let X be either the Euclidean space M n , the n- dimensional hyper- 
bolic space H n , or the unit sphere S n in R n+1 . We denote by 5 the set 
of all fc-dimensional totally geodesic sub manifolds of X, 1 < k < n — 1. 
The totally geodesic Radon transform of a function / on X is defined 
by 

0R/)(O = f f(x)d(X, (GH, (1.3) 
t 

where d^x stands for the corresponding canonical measure; see [H] for 
details. Let d(x,£) be the geodesic distance between x G X and (eH. 
We introduce the distance function 

( d(x,£) if X = W 1 , 
p(x,C) = < sind(x,0 ifX = § n , (1.4) 
( sinh d(ac, if X = W 1 , 

and the corresponding shifted dual Radon transform 



(R* r <p)(x)= J ^(0^(0. ^I, r>0. (1.5) 

The last integral averages <p(£) over all ^ G H satisfying p(x,£) = 
r. After pioneering papers by P. Funk and J. Radon, it is known 
[GGGl El IRu02al IRu02bl [Ru04] that / can be explicitly reconstructed 
from (p = Rf by applying differentiation of order k/2 (if k is odd this 
differentiation is fractional and interpreted in a suitable sense) to R*ip 
in the r 2 -variable. However, it was recently proved by Helgason [HI p. 
116] that in the case X = M. n and k even one can differentiate in r, 
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rather than in r 2 , and obtain the following result: 

f(x) = c k [d k r (R* r Rf)(x)] r=0 , (1.6) 

where Ck is a constant depending only on k. The constant Ck was not 
evaluated in [Hj. 

We compute this constant explicitly and generalize (jl.fip for arbitrary 
constant curvature space X and arbitrary 1 < k < n — 1. 

Notation. For the sake of simplicity, throughout the paper, we 

assume / to be infinitely differentiable. If X is noncompact, we also 

assume that / is rapidly decreasing together with derivatives of all 

orders. The letter c stands for a constant to be specified in every 

occasion; cr n -i = 2ir n ^ 2 jY{n/2) is the surface area of the unit sphere 
§n-i in R n 



2. Main Results 



Let p = p(x,£) be the distance function (ll.4p . 1 < < n — 1. For 
r > we denote 



(£»(*) = / ^(Op fc+1 " n sgn(p-r)^ 



(2.1) 



(L»(^) = / y>(0 P fc+1 " n log |p 2 - r 2 | de. 



(2.2) 



Theorem 2.1. Lei <p = i?/. 

(i) // is even, then 



d k r + \L*^){x) =d x f(x), 



r=0 



(2.3) 



where 



d 



x 



2 (-l)(fc+2)/2 (Tn _ fc _ i(Tfc _ i ( k _ ^ if X = R n ;H n ; 



2o- n _ k _ 1 a k a k „ 1 (k - l)\/a n if X = S n 
(ii) If k is odd, then 



d k +\L* rV )(x) =dxf{x), 



r=0 



(2.4) 



where 



( ' x " ^ 2n(-lY k - i y 2 a n ^k-iO- k a k - 1 (k-l)\/a n if X 
Theorem 2.2. Let (f = Rf . If k is even, then 

d k r \ x (r)(R* r v)(x) =c x f(x), 

r=0 



if X = R n ,U n , 



(2.5) 
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where 

( 1 if X = R n , 

\ x (r) = < (1 - r 2 )^- 1 )/ 2 if X = S n , 

I (l +r 2)(fc-l)/2 ^ X = U n , 

_ f {-lfl 2 (k- l)!<r fe _i »/ X = R", H", 
Cx -\2{-l) k '\k-l)\a k ^ ^/X = §^ 

Theorem 12.11 includes Mader's result for X = M. n , k = n— 1. Theorem 
12.21 gives precise form to the afore- mentioned result by S. Helgason and 
extends it to any constant curvature space X. 

The following lemma plays the key role in our consideration. 

Lemma 2.3. Let k be an odd positive integer, 

i 

^ k {u) = J (l - t,2)*/2-l ^ | u 2 _ w 2| ^ ( 2 _g) 


T/ien 

/ / \ _ / p k-i{u) if < u < 1, , s 

" { + ,(-l)(*-^9(u) ./ u > 1, [2 - 7} 

where Pk-i(u) is a polynomial of degree k—1, Q(u) = f^ [v 2 — l) k / 2 ~ l dv . 

This statement is presented in Mader's paper [M] in a slightly differ- 
ent form. We obtain (12. 7p as a consequence of a more general result, 
that might be of independent interest. 

Sections 3,4,5 contain the proof of Theorems 12.11 and 12 . 21 for the cases 
X = M n ,§ n , and EI™, respectively. Lemma [2.31 is proved in Section 6. 

3. The case X = W l 

We recall basic definitions. Let G n ^ (1 < k < n — 1) be the 
Grassmann manifold of fc-dimensional linear subspaces ( of W 1 and let 
5 = Q n ,k be the affine Grassmann manifold of all non-oriented fc-planes 
£ in M. n . Each /c-plane £ is parameterized by the pair ((,u), where 
C G Gn,fc arid it G (the orthogonal complement to ( in M n ). The 
manifold Q n ^ will be endowed with the product measure d£ = d(du, 
where d( is the SO (n)-invariant measure on G n ,u of total mass 1, and 
du denotes the usual volume element on £ . 

The A;-plane transform (i?/)(£) of a function / on X = W" - is defined 

by 

0R/)(O = Jf(u + v)dv, £=(C, (3.1) 

c 
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Throughout this section we assume that / is a rapidly decreasing in- 
finitely differentiable function. 

The shifted dual fc-plane transform is defined by 



(Kv)(x)= J ¥?(7Co+x+r 7 e n )d7, (3.2) 

SO(n) 

where (o is an arbitrary fixed dimensional subspace of M n , and e n = 
(0, . . . , 0, 1). This integral averages <p over all A;-planes at distance r 
from x. The case r = corresponds to the usual dual Radon transform 

m- 

The spherical mean of a locally integrable function / at a point 
x G M n is defined by 

(M t f)(x) = — [ f( x + te)de, t>o, (3.3) 



Cn-1 

Sn-1 



so that \mv(M t f){x) = f(x). 

Lemma 3.1. The following equalities hold: 

oo 

(R* r Rf)(x) = a^ x [ {M t f){x){t 2 - r^'Hdt; (3.4) 



^(Oo(p(x,0)de = ^-*-i / r n - k - l a(r)(R^)(x)dr. (3.5) 

Qn,k 

It is assumed that f, ip and a are measurable functions, for which either 
side of the corresponding equality exists in the Lebesgue sense. 

Equality (13. 4p can be found in |H, p. 115, formula (8)]. Both equal- 
ities were proved in [Ru04, formulas (5.2), (2.29)]. 
Proof of Theorem [2J] (the case X = M n ). We apply ([33} to 
(p = Rf and replace R*(p = R*Rf according to (13. 4p . Changing the 
order of integration, we obtain the following expression: 

oo 1 

o-n-k-io-k-i [ (M t f){x)t k dt [ {tvY-^aitv^l -v^-Hv. (3.6) 





If k is even, we choose a(-) in the form a(p) = p k+l ~ n sgn(p 2 — r 2 ). 
Then (jSSJ yields 

oo 

(L* r cp)(x) = (Tn-k-KTk-i [ {M t f){x)t k ^{r/t)dt, (3.7) 
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ip{u) = J sgn(v - u) (1 - v 2 ) k/2 ^dv. (3.8) 
o 

Setting 

1 u 
J (I _ V 2)fc/2-l du> e(u ) = y ( v 2 _ l)fc /2-l du> 



1 

Cfc = 



we have 



r- c , + 2(-i)'/>e M if o <«<i, 

rv 7 [ -c fc if u > 1. v 7 

This allows us to represent o^^o"^ (#) as 

r oo 

-c fc y (M t f)(x) t k dt + J (M t f)(x) t k [-c k + 2(-l) fe / 2 0(r/t)] eft 

r 

r 

= 2(-i)( fc + 2 )/ 2 /" (.M t /)(x) t k Q(r/t) dt 



oo 

+ y (M/)(x) t fc [-c fc + 2(-l) fc / 2 6(r/t)] eft. 
o 

Since J °° is a polynomial of degree k — 1, then 

r 

<9 r fc (L»(:r) = 2 ( x n _ fc _ l( 7 fc _ 1 (-l)( fe+2 )/ 2 ^ y (W)(x)t fe e(r/t)dt 

o 

r 

= 2 ( x n _ fc _ 1( 7 fc _ 1 (-l)( fe+2 )/ 2 ^- 1 y (M t f)(x) t k [d r Q{r/t)\ dt 



= 2a n ^ 1 a k ^{-l) {k+2 ^ 2 d k -\Kf){x), 

r 

(A r f)(x) = J (M t f)(x) (r 2 -t 2 f 2 -Hdt. (3.10) 
o 

We write the last integral as A 1 + A 2 , where 

r 

Ai = /(x) J(r 2 -t 2 ) k / 2 -Hdt = jf(x), 
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A 2 = (r 2 - t 2 ) k ' 2 - l [{M t f){x) - f(x)]tdt = r k h(r), 



h{r) = I (1 -t 2 fl 2 - l [{Mrtf){x) - f(x)]tdt. 

This gives 



kW(A r /)(>) = (k -l)\f(x), (3.11) 
i — >o 



and therefore, 



]hnd k+1 (L* r cp)(x) = 2{-lf + ^ 2 a n ^ l a k ^{k- l)\f(x), (3.12) 

r— w 

as desired. 

If k is odd, we write (13. 6 p with a(p) = p k+1 ~ n log |p 2 — r 2 | in the form 
Gri—k—i&k—i {A + -B(r)), where 

oo i 

* = **** J w- l *>- 



B(r) = J (M t f)(x)t k Mr/t)dt, (3.13) 
o 

and ipk being a function (12. 6p . By Lemma [2. 3 j 

r 

B{r) = [ (M t f)(x) t k [P k -!(r/t) + Tr(-l)^ 2 e(r/t)\ dt 



o 

oo 

k 



- I (MJ)(x)t k P k ^(r/t)dt 

r (M t f)(x) t k P fc _ x (r/t) dt+ 7 r(-l)( fc - 1 )/ 2 l{M t f)(x) t k Q(r/t)} dt. 



r 

oo 



The integral J °° is a polynomial of r of degree k — 1. Denoting 

(L^)(x)= J V (Op k+1 - n \og\p 2 ~r 2 \d^ p = p(x,0, (3.14) 



Gn.fc 
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and differentiating k + 1 times, we obtain 

r 

= tt(-1 ) {k - 1)/2 a M a k ^ d k r J (MJ)(x) (r 2 - t^'Hdt. 

o 

The last integral is already known; see 1 13. 10p . Hence, by 1 13. lip . 
lim9 r fe+1 (Z»(x) = cf(x), c = A-l^-^a^^a^ (k - 1)!, 

r— >0 

which completes the proof. 

Proof of Theorem 12.21 (the case X = W 1 ). For ip = Rf, equal- 
ity dH yields = afc-i [(-l) fe / 2 (A r /)(x) + (A r /)(x)], where 
(Arf)(x) is a function 113.101) and 

oo 

(A r /)(x) = J{M t f){x){t 2 -r 2 ) k ' 2 -Hdt. 

Q 

Since (A r f)(x) is a polynomial in the r- variable of degree k — 2, owing 
to 113.111) . we get 

hm^(i? T V)(x) = (—l) k/2 (k — 1)! cr fe _i f{x). 
i — >o 

4. The case X = § n 

We recall that § n is the unit sphere in M" +1 , d(-, •) denotes the ge- 
odesic distance on S n , and H stands for the set of all fc-dimensional 
totally geodesic sub manifolds £ of S n . Each £ G 5 is an intersection of 
S n with the relevant (k + l)-plane through the origin. Thus S can be 
identified with the Grassmann manifold G n+ i^ + i. The totally geodesic 
Radon transform of a function / on S n is defined by 

Wm = J f{x)dtx, £G ~, (4.1) 
t 

where d^x stands for the usual Lebesgue measure on £. Throughout 
this section we assume that / is an even C°° function. 

We set R n+1 = R k+1 © W n ~ k , R k+1 = Re 1 © . . . © Rejt+i, 
M. n ~ k = © ... © Ke n+ i, e$ being the coordinate unit vectors. 

Let £ = S fc be the unit sphere in R fe+1 . Given x G § n , £ G S, 
we denote by an arbitrary rotation satisfying r x e n+ \ = x and set 

Let .K" = SO(n) be the group of rotations about the x n+ \ axis. 
For G [0,7r/2], let gg be the rotation in the plane (efe + i, e n+1 ) with 
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the matrix 



sin 9 cos 9 
— cos 9 sin 9 
function <p on H is defined by 



. The shifted dual Radon transform of a 



(i?»(x)= / <p x (P9e%)dp. (4.2) 



A" 



The case = corresponds to the usual dual Radon transform [H] . 
A geometric meaning of (i?gV 9 )( x ) is as follows. Let 

xe = gee n +i = e k+ i cos 9 + e n+1 sin 9 

so that d(xo,£o) = 9. Hence, 

(R* e <p)(x)= J <p(Odm= J <p(0dm = (R»(x)i (4-3) 

where r = sm9, p(x,£) = smd(x,£,), dp{£) is the relevant normalized 
measure; cf. (II. 5ft . 

We need one more averaging operator. Given x G § n and t G (—1,1), 
denote 

(I _ t 2)(l-n)/2 r 

(M t f)(x)= { - 1 / f(y)da(y). (4.4) 

{y& n : x-y=t} 

The integral (I4.4p is the mean value of / on the planar section of § n by 
the hyperplane x ■ y = t, and da(y) stands for the induced Lebesgue 

measure on this section. Thus, lim(M t /)(x) = fix). 

t— »i 

Lemma 4.1. Let p(x,£) = sin[d(:r, £)]. The following equalities hold: 



{RlRf){x) = 2a k _ 1 [ 
Jo 



l-r 2 ) k / 2 -\M TCOse f)(x)dr; (4.5) 



a(pMMZ)d£ (4.6) 



i 

0~n-k-lO~k 





1 _ r 2 )( fc - 1 )/ 2 r n - fc - 1 a(r)( J R* n - lr ^)(x) dr. 



It is assumed that f, ip and a are measurable functions, for which either 
side of the corresponding equality exists in the Lebesgue sense. 

These equalities were proved in |Ru02bt pp. 485, 480]. 
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Proof of Theorem EH (the case X = § n ). Let <p = Rf. We write 
( 14. 5 p in the form 

(R* e <p)(x) = ?^=i \\p 2 - s 2 ) k ^-\M s f)(x) ds, p = cos 6 £ (0, 1). 
P Jo 



Setting r = sin 9 = y/l — p 2 , s = y/1 — t 2 , and denoting 

(MJ)(a:) = (l-t 2 )- 1 /2(M vCT /)(x), 
so that lim(M(/)(x) = f(x), we obtain 



t-+0 

•1 



(RW)(x) = n 2 ?JLiw9 / (M t /)(x)(t 2 -rY /2 -^t. (4.7) 



(1 _ r 2)(fc-l)/2 

Then we apply ( 14. 6 p to ip = i?/ and replace R* ^ <p according to (14. 7p . 
Changing the order of integration, we obtain 

J a{p{x,0)cp(0dt (4.8) 



(M t f)(x)t k dt / {tvf-^aitv)^ - v 2 ) k,2 - l dv, 



i 

= c 



c = . 

This mimics (I3.6P and we continue as in Section [3j 

If is even we choose a(-) in the form a(r/) = 77 fc+1_n sgn(r] 2 — r 2 ). 
Then for <p = Rf, (USD yields 

i 

= c J (M t f)(x)t k ij(r/t)dt, (4.9) 

Q 

where ip is a function from the previous section; see ( 13.81) . (13.91) . Hence, 
as in Section 3, (L*<p)(x)/c can be written as 

r 1 

-c k J(M t f)(x)t k dt + J(M t f)(x)t k [-c k + 2(-l) k / 2 Q(r/t)]dt 

r 



2 (_i)(fe+2)/a / (M t f)(x)1*B(r/t)dt 



o 

i 



+ / (MJ)(x)t fc [-c fc + 2(-l) fc / 2 6(r/t)] ( it. 



o 
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Since J (. . .) is a polynomial of degree k — 1, then verbatim of the 
corresponding reasoning from Section 3 yields 

\imd k+ \L* rV )(x) =c(k- 1)! f(x) = 2 °n-k-i°k°k-i (Jfe _ 1} , 

r-+0 a n 

If is odd, we choose a(-) in the form a(r/) = rj h+1 ~ n log |r/ 2 — r 2 | and 
write (14 .8p as c(A + B(r)), where 

i i 
A = 2 J (M t f)(x)t k hgtdt J {I- v 2 ) k/2 ~ l dv = const, 



o 

i 



B(r) = J(M t f)(x) t k Mr/t) dt, 
o 

and being a function (j2.6P ; cf. (I3.13p . By Lemma [2~3| 

r 

B(r) = [ (M t f)(x) t k [P k -x{r/t) + n(-l)^ 2 e{r/t)] dt 



o 

i 

k 



+ (M t f)(x)t k P k ^(r/t)dt 



r 

1 



t /)(z) t fc P fc _i(r/t) dt+ 7 r(-l)( fe - 1 )/ 2 / (M t /)(x) t fc 9(r/t) 



The integral Jg 1 is a polynomial of r of degree k — 1. Denoting 

(3V)(aO= /^(Op fc+1 " n log|p 2 -r 2 |^, P = (4-10) 



(cf. (I3.14p ) and differentiating k + 1 times, we obtain 

r 

d r fe+1 (Z»(x) = cvr(-l)( fc - 1 )/ 2 9 r fc / (M t /)(x) (r 2 -t 2 ) k/2 -Hdt. 



o 

As in Section [3], this gives 

\imd k+1 (L* r <p)(x) = cf(x), (4.11) 

2vr (_l)(fc-i)/2 an _ fc _ 1 a fc a fc _ 1 (A;-l)! ~ 
c = (= d x )- 
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Proof of Theorem 12.21 (the case X = E> n ). Let ip = Rf and write 
( 14. 7p in the form 

where 

1 r 

h(r)=J(M t f)(x)(t 2 -r 2 ) k / 2 -Hdt, I 2 (r)=J(M t f)(x)(r 2 -t 2 ) k / 2 -Hdt, 
o o 

I\{r) being a polynomial of degree k — 2. Hence, as in (13.111) . we obtain 

lim#[(l - r 2 )^' 2 (R^Mx)] = 2(-l) fc /V fc _ 1 (k - 1)! f(x). 

i — >u 

5. The case X = U n 

Let E 11 ' 1 , n > 2, be the pseudo-Euclidean space of points x = 
i, . . . , x n+ i) G M n+1 with the inner product 

[x, y] = -xiyi - ... - x n y n + x n+1 y n+1 . (5.1) 



[Xl 



We realize the n-dimensional hyperbolic space X = H" as the upper 
sheet of the two-sheeted hyperboloid 

W 1 = {xe W 1 ' 1 : [x, x] = 1, x n+1 > 0}. 

Let H be the set of all k- dimensional totally geodesic submanifolds 
£ C H n , 1 < k < n — 1. As usual, ei, . . . , e n+ i denote the coordinate 
unit vectors. We set JJT+ 1 = R n ~ k © R fc+1 , where 

R n ~ k = Rei © . . . © Re„_ fc , R k+1 = Re r ^ k+1 © ... © Re n+ i, 

and identify M. k+1 with the pseudo- Euclidean space E k ' 1 . 

In the following x = (0, . . . , 0, 1) and £ = H n n E M = H fe denote 
the origins in X and H respectively; G = SO (n, 1) is the identity com- 
ponent of the pseudo-orthogonal group 0(n, 1) preserving the bilinear 
form (15.11) ; if = SO{n) and If = SO(n—k) x SO (k, 1) are the isotropy 
subgroups of x and £ , so that X = G/X, H = G/H. One can write 
/(x) = f(gK), (p(£) = <f(gH), g E G. The geodesic distance between 
points x and y in X is defined by d(x,y) = cosh _1 [x, y}. 

Given x G X, £ G S, we denote by r x ,rg(G G) arbitrary pseudo- 
rotations such that r^-Xo = x, r^ = £, and write 
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The totally geodesic Radon transform (i?/)(£) of a rapidly decreasing 
C°° function / on EI n is defined by 



W)(0= / f(x)d& 



f^x )dj, 



(5.2) 



SO (fc,l) 



As is [BR], let 



9e 



cosh 9 sinh 9 

I n -1 

sinh 9 cosh # 



where i n _i is the unit matrix of dimension n — 1. The shifted dual 
Radon transform of a function p on 5 is defined by 



(5.3) 



K 



The case 9 = corresponds to the usual dual Radon transform. 

A geometric meaning of {R* d p){x) is as follows. Let d(x, £) be the 
geodesic distance between i6l and £ G S. We set 

x<9 = (76»e n +i = ei sinh # + e n+ i cosh 6* 

so that d(x0,£o) — d{xe,e n+ \) = 9. Then 

d(x,r x -fg e %) = 9 (5.4) 

for all 7 G if and x G X. Thus, 



(£»(:r) 



v(0 MO 



<p(Z)dn(Z) = (ISip)(x), (5.5) 



p(a:,£)=r 



where r = sinh#, p(x,£) = sinh d(x,£), dp,({;) is the relevant normal- 
ized measure; cf. (11. 5ft . If ip G -^ 1 (2), then 



oo 

y ^(0 rf£ = <x n - fe -i y (R* e <p)(x) dv{9) ^x G A, 

3 

dv{9) = {sinh 9) n - k - 1 {cosh 9) k d9; 
see (EEl formula (2.30)]. Given i G I" and t> 1, let 

(t 2 - l)^-™)/ 2 



(5.6) 



/(y)d<%) (5.7) 



be the spherical mean of / on A = HP, where da{y) stands for the 
induced Lebesgue measure. Then lim(M t /)(x) = f{x). 

t-¥l 
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Lemma 5.1. Let p(x,£) = sinh [d(x, £)]. The following equalities hold: 

oo 

(RlRf)(x) = ^± j{M s f){x){s 2 -r 2 ) k / 2 - l ds, r = cosh£; (5.8) 



a(p(x,0)<p(Od£ (5.9) 



0~n 



n-k-i / (l + rT-^r^airXR* ,^p)(x)dr. 



o 



It is assumed that f, (p and a are measurable functions, for which either 
side of the corresponding equality exists in the Lebesgue sense. 

Proof. Equality (15.81) can be found in [H]; see also |Ru02aj formula 
(2.3)]. Equality Q follows from (ESD, (E2D, and Indeed, 



a(p(x, 0) f(0 d£ = I 0(0 d£, 0(0 = a(p(e n+h 0) f x (0- 

Hence, 

I = <T n -k-i I (R*e0)(x) du(9) 





oo 



o-n-k-i / dv(6) I a(sinhd(e n+1 ,^g 6 1 ^ ))<p x (jg d 1 ^ ))d"f 



K 

oo 



= cFn-k-i J a(smh.6)(R* e (f)(x)dv(6). 



This gives (ES). □ 

Proof of Theorem O (the case X = H n ). Let 93 = Rf. Setting 
r = sinh 9 = y/r 2 — 1, s = \/I + 1 2 , and denoting 

(M t f)(x) = (1 + i 2 )- 1 ^^/)^), (5.10) 

so that \im(M t f)(x) = f(x), we obtain 



r-OO 
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Then we apply ( 15. 9 p to <p = Rf and replace R*. mh -i <f according to 
( 15. lip . Changing the order of integration, we obtain 

a(p(s,0M0# (5-12) 

oo 1 

= a n - k -i* k -if(M t f)(x) t k dt j (fe)"- fc " 1 a(^)(l -v 2 ) k / 2 - l dv. 



This mimics H3 .61) and we proceed as in the previous sections. Specif- 
ically, choose a(-) in the form a(rj) = r] k+1 ~ n sgn(rj 2 — r 2 ). Then for 

<p = Rf, (JS32D yields 

oo 

(L* r ip)(x) = o rM o k _ x J (M t f)(x)t k ^(r/t)dt; (5.13) 

o 

cf. (14. 9p . Hence, as in Section 3, 

limd r fc+1 (L^)(x) = 2(-l)( fc+2 )/V n _ fc _ 1 a fc _ 1 (k - 1)! f(x). (5.14) 

If k is odd, we choose a(rj) = r] k+1 ~ n log \rf — r 2 \ and set 

= j ^)p k+1 ' n \og\ P 2 -r 2 \d^ p = p(x,0, (5-15) 

as in (14. lOH . This gives 

lim9 r w (L»(x)=c/(x), c = 7r(-l)( fe - 1 )/ 2 a n _ fc _ 1 a fc _ 1 (A;-l)!. 

i — 

Proof of Theorem 12.21 (the case X = HP). For ip = Rf, equality 
(Oil yields 

(^W^X*) = (1+ ^)"(Li)/2 [(-l) fc/2 (AJ)(x) + (A r /)(x)], (5.16) 
where (cf. (I5TT0|) . (|3TT0|l ) 

r 

(A r f)(x) = [ (M t f)(x) (r 2 - t 2 ) k l 2 ~Hdt, 



o 



(A r /)(x)= {M t f){x){t 2 -r 2 ) k / 2 -Hdt. 



Hence, as in ( 13. lip , we obtain 

lim^[(l + r 2 )^' 2 {Rl^ rV ){x)\ = (-l) k/2 a k -i (k - 1)! f(x). 
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6. Proof of Lemma [2731 

The statement will be obtained as a consequence of the following 
more general result, which generalizes the reasoning from [M] . 

Lemma 6.1. Let m G Z, — 1 < a < m + 1, a ^ < u < oo. Then 

i 

: j(l+O a (l-O m - a log\Z-u\dt (6.1) 



1* 



where 



li a (u)Q a (u) + P m+1 (u), (6.2) 

-7T cot«7r, if < u < 1, 

— l) m 7r csc an, if 1 < u < oo, 

u 

Q a (u) = J(l + ^\1-^-^, (6.3) 
l 

and P m+ i(u) is a polynomial of degree m + 1 defined by 

m+l 

P m+ i(u) = 0(1) - (-l) m 7r cscavr ^ A r (u r - 1), (6.4) 



r=l 

m+l— r 



m — a\ a 



r / — ' \ £ / \m + 1 — r 

^=i v / \ 1 

Proof. The result will be obtained using the classical tools of complex 
analysis. Consider auxiliary functions 

C(z) = (z + l) a (z - l) m ~ a , h{z, u) = C(z) log -; z G C. 

z — 1 

To fix a single branch of C(- 2 )) we cu t the z-plane along the segment 
— 1 < < 1 and choose arg(z ± 1) G [0, 2n]. Then ((2) ~ z m , z — >■ oo, 
and 

r (-i)™ e Ti "|l + C| a |l-^| m - a , if-l<£<l, 
C(£±*0)=<^ (6.5) 
[ |1 + £| Q |1 -^| m " a , if 1 < f < oo. 

Case 1. Let 1 < u < oo. Continue the cut [—1,1] along the real 
axis up to the point u and fix a branch of the logarithmic function so 
that arg(;z — u) G [0, 2tt] and, as before, arg(z — 1) G [0, 2ir]. Then 

z — u 

log ~ 0, z — > oo, 

2 — 1 
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and 



log = log 

6 f ± iO - 1 s 



27T, 1 < £ < W, 

o, o<e<i. 



(6.6) 



At the point 2; = 00, the function h(z,u) has an integer order, and its 
residue is 



c_i = res h(z,u) = — : / h(z,u)dz. 
z=oo 2ni J 



(6.7) 



G- 



Here c_i is the coefficient at z~ x in the Laurent series of h(z,u) in 
a neighborhood of the infinite point, and C~ is a circle of radius R 
centered at z = 0. It is assumed that i? is big enough, so that z = 00 
is the only singular point outside of C~ . The direction of integration 
on C~ is chosen clockwise. Deforming the contour C~ , we obtain 

res /i(£, u) 



2vri 



/i(f,tt) d£, 



(6.8) 



where L = [u, 1]- U [1, -1]" U [-1, 1]+ U [1, u] + , and [a, b] + and [a, b]~ 
are the upper and lower sides of a cut [a, b], respectively. Owing to 
(ESD and fl62D, (E3D yields 

• 1 



res h(z, u) = Q a (u)- 



sin air 



71 



;i+o Q (i-o m_a iog 



d£. 



In other words, 

0(«) = 0(1) + /^(m) © a (w) — ^ a (u) res h(z,u), (6.9) 

z=oo 

and it remains to evaluate the residue on the right-hand side. We set 
h(z,u) = -W(z,u) + W(z,l), (6.10) 

where 



W(z,u) 



—z 



1 + 



<h z m ~ p 

p=0 



1 - 



EVE 

q=0 r=l 

Vi — Q 
9 



log 1 



r 



A simple calculation gives W(z,u) = 72 1 + terms of degree 7^ — 1, 
with 



ro+1 



;E (-D'( m 



7 = $^A r 



A, 



r=l 



a 

m + 1 — r 
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Hence, 
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m+1 

res h(z,u) = V] A r (w r - 1) = Q m +\{u) 

z=oo ' ' 

r=l 



(6.11) 



and, by ( 16. 9p . we finally get 0(m) = fi a (u) Q a (u) + P m+ i(u), where 



P m +i(u) = 0(1) - fi a (u) Q m+ i(u) = 0(1) 



-l) m 7T 



sm air 



Q m+l {u). (6.12) 



Case 2. Let < u < 1. Then the substitute for (I6.8P has the form 

res h(z,u) = ^— [ h(£,u)d£, (6.13) 

z=oo ZTIZ J 

L' 

where V = [l,u]~ U [u, — 1]~ U [— l,u] + U [u, 1] + , and the residue of 
h(z,u) at the infinite point remains unchanged. Now, (16. 3p . (16.51) . the 
counterpart of formula (16. 6p for this case 



i , 



£ ± iO - 1 
and ([SIT]) yield 



^-1 



± 



?7T, U < £ < 1, 

0, < £ < u, 



(6.14) 



<2m+l(w) 



•I)™ - 1 COS «7T 6 a (w 

-l) m sina7r f 1 



71 



(i + Q (i-O m "°iog 



e-1 



This gives 



(«) = 0(1) — 7T COt COT Q a (u) - 

= /i a (u)Q a (u) + P m+ i (u) 



-l) m 7T 



sm an 



Qm+i(u) 



as desired. 



□ 



Let us turn to the particular case a = m — a = | — 1, where A; is an 



odd positive integer. Then 



y (1 - ^ a )*/2-i ^ | f _ „[ df, (6.15) 



-1 



and Lemma [6.11 yields 

6(u) = P k -i(u) + 



0, if0<«<l, 

7r (_l)(fc-l)/2 ( u ) ) if I < u < QO) 



(6.16) 
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where Pk-i{u) is a polynomial of degree k — 1, and 

u 

e(n) = J{e-i) k/2 - l dt 

1 

Since 

1 1 

= (/+/)(•••) = /a - e 2 )^ 2 - 1 log ie 2 - « 2 i 

-10 

we arrive at Lemma [2.31 
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